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Abstract 

Most of theoretical physics is based on the mathematics of func- 
tions of a real or a complex variable; yet we frequently are drawn to try 
extending our reach to include quaternions. The non-commutativity 
of the quaternion algebra poses obstacles for the usual manipulations 
of calculus; but we show in this paper how many of those obstacles 
can be overcome. The surprising result is that the first order term in 
the expansion of F{x + 5) is a compact formula involving both F'[x) 
and [F{x) — F[x*y\/{x — x*). This advance in the differential calculus 
for quaternionic variables also leads us to some progress in studying 
integration. 
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1 Introduction 



We are very familiar with functions of a real or complex variable x which we 
can expand, in the mode of differential calculus, as 

F{x + 8) = F{x) + F'{x)5 + ]^F"{x)5^ + (1.1) 

But what if we consider a quaternionic variable 

X = xo + ixi + ix2 + fcxs, (1.2) 

involving four real variables, x^, /i = 0, 1, 2, 3, along with those quaternions 
j, k which do not commute with one another. 

= = k"^ = —1, ij = —ji = k, etc. (1.3) 

The small quantity 6 will also involve all those quaternions. How then can 
we expect anything as neatly packaged as Equation fll.ip ? 

This is a long-standing challenge to mathematicians. The earliest at- 
tempt to extend the usual concept of the derivative, with a quaternionic 
dx, failed dramatically. The subsequent approach focused on the four real 
variables, 

dFix) = E ^dx,. (1.4) 

That approach, often called quaternionic analyticity, springs from the work 
in the 1930's by R. Fueter [1] and his school, with more accessible articles 
reviewing that subject available in references [2], [3]. Some more recent 
attempts to advance that work may be found in references [1], [S], [H]- In 
Appendix C I provide a rough summary of the Fueter approach. 

The first new result presented in this paper is an alternative approach 
to the differential calculus - something between relying on the whole quater- 
nionic variable, and resorting to the four-component real variables, as 
in Eq. fll.4p . This starts, in Section 2, with the separation of the quater- 
nionic displacement 6 into two parts, one "parallel" and the other "per- 
pendicular" to the quaternionic variable x as it may be envisioned in that 
four- dimensional space. 

The subsequent Sections show how this leads to a surprisingly compact 
and general formula for the quaternionic version of the expansion fll.ip : 

F{x + 6)= F{x) + F\x) 5\\ + (F(s) - F{x*))/{x - x*) 5^ + 0{5^). (1.5) 
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In the second part of this paper we look at integration; and find that the 
new form of the quaternionic differential leads to new results in this other 
realm of calculus. 

2 Local Coordinates 

The standard approach to quaternionic variables starts with a global set of 
imaginary coordinates; 

X = xo + ixi + jx2 + kxs, 5 = 60 + i6i + j62 + kd^. (2.1) 

We now want to write x in a different way: 

x = Xo + Uxr; r = xj + xj + x^, ul = -1 (2.2) 

where Ux is a unit imaginary that varies in its i,j, k composition as x moves 
from one point to another in that 4-dimensional space. This is analogous to 
the use of polar coordinates in 2-dimensional Euclidean space. 

Now we want to decompose the quaternionic quantity 5 in a particular 
way that refers to this local coordinate system. 

S = S\\+6±, 6\i = ^{6 -Ux 6 Ux), S± = ^{6 + Ux 6 Ux), (2.3) 

which leads to the algebraic relations 

S\\ Ux = Ux S\\, S± Ux = —Ux S±. (2.4) 

The essence of this approach is expressed in the nomenclatures "parallel" 
and "perpendicular" for these two components of 6 as they relate to the 
local quaternion x. The most useful way to write these relations is 

6\\ X = X 6\\, 6± X = X* 6±, (2.5) 

where * is the complex conjugation operator that changes the sign of all 
imaginaries. 

Now we shall give three examples of how to expand F{x + 6) with this 
simple machinery. 

^Thc efBcacy of this technique was discovered as the resuh of a more long-winded 
calculation, which may be seen in reference [7]. 
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3 The Function F{x) = x 

We calculate directly, 



n 



n-l 



{x + 5)" = + ^ x"-"-^ 6x"' + 0{6^). (3.1) 



m=0 



Putting in the separation 6 = 6\\ + 6± and using the properties of Eq. fl2.5p . 
the sum becomes 



n-l 



^ = ^ (x^-i S\\ + x"-"^-^ X* "'Si_)=n x"-^ S\i + (x" - X* ")(x - x*)-^ S^, 

(3.2) 

where we evaluated a finite geometric series. 

4 The Exponential Function 

For a general quaternionic variable x, we can define the exponential function 
in the usual way: 

e^= lim(l + ^)^ (4.1) 
and this leads us to the expansion, 

g(x+5) = e"[l + f' ds e-"" 5 e"" + 0{5% (4.2) 

which is correct for the situation where x and 5 do not commute. For a 
derivation of this formula, see Appendix A. 

Putting in a real parameter p, and following the course set above, we get 
the expansion, 

^p(x+&) _ = dsp e^'~'^P%5ii + 5x)e'^" = (4.3) 

ds p e^^ + dsp e(^-")P^ e'^^* 5± = (4.4) 
Jo Jo 

p e^" 5|| + (e^" - eP"*)(x - x*)-^ 5± (4.5) 

to first order in 6. 
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5 General Analytic Function F(x) 

For a general analytic function F{x) of a quaternionic variable x, we start 
by assuming a representation as a Laplace transform: 

F{x) = Jdp f{p) e^^ (5.1) 

where p is a real variable. We then use the result of the previous Section to 
obtain 

F{x + 5)- F{x) = F\x) 5\\ + {F{x) - F{x*)) {x - x*)-^ 6^ + 0{6^), (5.2) 

where F'{x) is the derivative of the function F{x) calculated as if x were a 
real variable. 

This is our general result. The particular result of Section 3, for F{x) = 
x"-, also fits this general formula; and thus it also works for any power series 

Fix) = En Cn X". 

The authors of reference [1] have taken an approach somewhat similar 
to what is done here. They introduced a local unit imaginary (which they 
call iota) that is the same as what we have defined as Ux- However, they 
limit their differentiations to displacements that are restricted to the two- 
dimensional space of what we call 6\\ without allowing any of 6±. In that way 
they merely reproduce what is known about ordinary complex variables. 

6 Further Exercises 

Let us define the first-order differential operator V, from Eq. fl5.2p . as 

F{x + 6) = F{x)+V F{x) + 0{6^) (6.1) 

with 

V F{x) = F'{x) 5\\ + {F{x) - F{x*)) {x - x*)"^ 5±. (6.2) 

Several interesting exercises are now suggested. 

Calculate V {F{x)G{x)) and verify the applicability of Leibnitz' rule. 
Calculate V and also V F{G{x)) . 
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7 Alternative Arrangements 



Still another way to represent our result for the first-order differential is in 
terms of some partial derivatives, defined as follows. 

^ , dF(x) , dF(x) , 

V F{x) = — ^ dx\\ + dx^, (7.1) 

OX\\ ox± 

^ = F'ix), ^ ^ (Fix) - Fix*))ix - XT'. (7.2) 

Suppose we restrict the functions F{x) to be real: that is, the coefficients 
Cn in F = J2n Cn x"' OT the amplitudes f{p) in the Laplace transform are real 
numbers. Then it is noted that the terms in Eq.( l5.2p can be written with 
the displacement quaternions, 6\\ and 6±, written either to the right or to 
the left of their accompanying factors. This is obvious in the case of F'[x), 
since 6\\ commutes with x. For the second term, we know that 6± does not 
commute with x; it takes the complex conjugate. But we note that the whole 
expression(F — F*) /{x — x*) is real; therefore this rearrangement is possible. 

The same rearrangement can be done with Eq. fl7.ip . 

These considerations lead us to note that the second term in the equation 
for T>F{x) can be written in terms of commutators as 

{F{x) - F{x*)){x - X*)-' = [C, F{x)], (7.3) 

where C is defined by 

[C,x]=5^, C = ^^S^. (7.4) 

X* — X 

What is somewhat surprising about this alternative arrangement, Eq. fl7.3p . 
is that the expression on the left hand side is manifestly non-local, involving 
things evaluated at the point x and also at the remote point x*; but the right 
hand side appears to be local, involving only x. This confusion is removed 
when one recognizes that C is a non-local operator, involving 6±, which 
changes x to x*. 



8 Second Order Terms 

Let's return to the exponential function fl4.ip and proceed with the expansion, 

fl rl rl-t 



e^[l+ / ds e-'''6e'''+ f dt f c/s e-(^+*)^(5e*^5e^^ + 0(53)]. (8.1) 
Jo Jo Jo 
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The best approach is to combine the exponential function and the Laplace 
transform from the beginnine:. Writing F{x + 6) = F{x) + F^^) + F^^) + . . 
we now look at 

Again, we decompose 5 and after a bit more work arrive at the result for the 
second order term, 

= 5j + {F{x) - F{x*)) {x - x*)-^ {6^6\\ - 66±) + 

F\x) {x - x*)-^ + F\x*) {x* - x)-^ 5^5\\. (8.3) 

It is also true that, with the first order term given as F^'^^x) = T> F(x), 
the second order result can be written as 

F^^\x) = ^VV F{x). (8.4) 

To verify this one needs the preliminary formulas, 

Vx = 5, V X* = 5*, V5 = 0, Vu^ = - 5±, (8.5) 

r 

along with 6u 6± = 6± 6\\ and 6]_ = —S±. See further in Appendix D. 



9 Discussion on Differentials 

It is surprising how simple and how general the new results obtained here 
are. It is also noteworthy that our differential operators are no longer local: 
they involve F{x*) along with F{x). 

One may ask what restrictions there are on the functions F{x) considered 
above. At first, one would say that they should be real analytic functions; 
having terms like xax where a is a general quaternion would certainly cause 
trouble]^ One can extend this condition slightly by allowing F(x) (but not the 
function G{x) in Section 6) to be a real function with arbitrary quaternions 
multiplying from the left. That is, the power series form F = J2n Cn x"' 
could have arbitrary numbers c„. 

^This use of the term "real analytic" differs from that found in reference [3]. 
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This bias to the left-hand side can be reversed if we change the original 
steps (14.21) . setting s 1 — s, and fl5.ll) . putting f{p) on the right-hand side. 

The Taylor series we have discussed above are expansions about the origin 
X = 0. In the usual complex analysis such power series may be about any 
fixed point x = x/; but such a quaternion constant put in the middle of our 
expressions would appear to cause trouble. That trouble could be avoided 
by limiting to be real; but there is a better way. If we define a new 
quaternionic variable y = x — Xf then we may proceed as done above only 
using the appropriate unit imaginary Uy, instead of the original Ux in order to 
separate the displacement S into "parallel" and "perpendicular" components. 

One may also ask if this general method may be applied to some other 
kind of non-commuting algebra beyond the quaternions. I believe that some- 
thing very similar can be done starting with a Clifford algebra. Other exam- 
ples are given in Appendix B and in reference [8]. 

10 Introduction to Integration 

When the conventional approach to analyticity of quaternionic functions 
failed in the differential calculus, the main push was then in the realm of 
integral calculus. 

The key result of the Fueter school was a third order differential equation 
that could define a "regular" function of a quaternionic variable, just as 
the Cauchy-Riemann equation was a first order constraint on functions of a 
complex variable z = x+iy. That approach is described roughly in Appendix 
C. Their result is a focus on integrals over a 3-dimensional surface in the 4- 
dimensional space. 

With the constructon of the quaternionic differential operator V we can 
do something quite different about integration, as is shown in the following 
two Sections. 

11 The Line Integral 

In ordinary calculus of functions of the real variable t, we know what is meant 
by an integral, such as / f{t) dt. But when we first consider quaternionic 
(or other non-commuting) variables it is unclear even how to write such an 
expression. We shall pursue that path in Section 12. 
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Alternatively, we can start with the defining relation between the integral 
and the differential: 

dfit)= f^dt = f{h)~f{a)- (11.1) 
Ja at 

and this is what we shall generalize for our non- commuting quaternionic 
variable 

r VF{x) = F{x,)-F{xa). (11.2) 

J a 

We define this integral as an additive operation along a path in that 
four-dimensional space of the real variables x^, 

where s is a real continuous parameter. 

Next, we subdivide that path, whatever it may be, into a large number 
of infinitesimal increments. 

= E / ' / = / ' n = l,...,N (11.4) 

where Xq = Xa and a; at = a^fe- 

In any segment of this path we choose the line of integration, with the 
integrand VF{x), to be the sum of two infinitesimal parts: 

Xn - Xn-l = S = 6\\ + 6±. (11-5) 

The first part is "parallel" to the direction of x at that point, giving the 
contribution 

^ VF{x) = F'(x)5||. (11.6) 
Then the second part is "perpendicular" , giving the contribution 

VF{x) = [F{x) - F{x*)]{x - x*)-^ 6^. (11.7) 

The sum of these two parts is thus nothing other than 

F{Xn)-F{Xn-l) (11.8) 

to first order in the interval 6. The entire sum then results in Eq. flll.2p . 
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Another general proof can proceed as follows. If we start with the coor- 
dinate along the path x{s) = Xpath{s), then we can simply write, 

Vx{s) = ds^^ (11.9) 
as 

since there is no commutativity problem in this representation. It is also true 
that we can express any function composed of powers of x as 

F{x{s)) = A{s) + B{s) x{s) (11.10) 

where A and B are real functions, the only quaternions being in the single 
factor x{s). We then see that the integral becomes quite ordinary: 

r VF{x{s)) = /' ds^^^^^ = F{x{s))\l = F{xb) - F{xa). (11.11) 

Ja Jo dS 

Since this differential operator V obeys the Leibnitz rule we get the iden- 
tity, usually called "integration by parts" , 

t F{x) VG{x) = F{x,)G{x,) - F{xa)G{xa) - t {VF{x)) G{x). (11.12) 

J a J a 

Loosly speaking, integration is the inverse of differentiation. What we see 
in Eqs. (Ill.ip and ( 111.2^ is one statement of that relationship. But there is 
also the other form, which is stated for real variables as 

I / m dt' = m. (11.13) 

For our quaternionic variables we start by looking at 

f V^. F{x') (11.14) 

and then apply the first differential operator to the coordinate x in two 
parts: first the 5\\ part and then the 5± part. The result is just the integrand 
evaluated at the point x\ 

= V^F{x)- (11.15) 

and this is just what we should expect from the right hand side of Eq. f lll.2p . 
with Xb replaced by x. 
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12 The Other Line Integral 

If we look at the common real integral and try to guess how to generalize it 
to the non-commutative quaternions, we might start with, 

J fit) dt ^ ^ J [dx F{x) + F{x) dx) ; (12.1) 

but why should dx only appear on the outside; why not also in the middle 
of the function F{x)l 

Let's try a most symmetrical arrangement with the function F{x) = x": 

Jrdt ^ I {dx x'^+x dx x'^-^+x'^ dx x''-^ + . . .+X'' dx). (12.2) 

But we can recognize that the long expression in parentheses on the right 
hand side of this is nothing other than Dx""*"^: 

T>F{x) = F{x + dx) — F{x), to first order in dx. (12.3) 

So we would then write, 

e dt / Pa;"+^ = , (12.4) 

n+lJ n+1 ^ ^ 

using our defining Eq. flll.2p . Now, this looks quite familiar. 

We can extend this to any power series and thus offer the following rule. 
For any analytic function of a real variable f{t), for which we know the 
integral, 

' f{t)dt = h{t), (12.5) 



we can make the correspondence to quaternionic integration as follows: 

J f{t)dt — > J Vh{x)=h{x). (12.6) 

While this may look trivial for real and complex variables, it is something new 
for non-commuting variables. This arises because we have carefully defined 
and studied the operator V. 
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13 Discussion on Integration 



Following what was stated earlier, we do require the functions F{x) to be real 
analytic functions along the path of integration. Terms such would 
be allowed only for real constants a. 

Our first new result Eq. flll.2p implies that the result of the integration 
depends only on the end points and is independent of the path. This is true 
if we also require that the function F{x) be single valued. Then, we have the 
result that the integral over any closed path, ending up at the same point 
where it started, is zero. This is a significant new result, carrying the world 
of contour integration over from the complex domain to the quaternionic. 

Our second new result, Eq's. fll2.5l) and fll2.6l) . opens up considerable 
possibilities for integration of quaternionic functions. 
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Appendix A - Expansion of the Exponential 

Here we give a derivation of the formula (14.21) for any non-commuting 
quantities x and 6. 

N-l r 

In taking the limit —>■ oo, we convert the sum over m to an integral over 
s = ^ and this yields 



e^' + ds e(i-^)^ 6 e^^ + 0{6^). (A.3) 
Jo 



Appendix B - SU(2,C) 

Here we shall extend the general method used above for a quaternionic 
variable to something built on a Lie Algebra - specifically SU(2). 
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Here is the Lie algebra: 

[Jl,J2] = Jz, [J2,h] = Jl, [J^,Jl] = J2. (B.l) 

where the three J's are understood to be matrices over the complex numbers. 
In particular we shall use the relations 

e^Ji j^e-^-^' = Jicos9+J2sin9, e^-^' J2 e~^-^'' = J2 cosO-Ji sinO, (B.2) 

which follow from (IB.ip . 

The new variable x is to be constructed with four real parameters as 

X = xo I + xi Ji + X2 J2 + xs J3 (B.3) 

and we want to expand F{x + 5) = F{x) + F^-^^ + 0(5^), where 5 is a small 
quantity in that same space of matrices as x. Our first step is to define a local 
coordinate system at the given point x. By a suitable linear transformation 
(rotation) of the Lie algebra we make the coordinate x appear as 

X = xo I + rJ-i (B.4) 

where we recognize that r"^ = xf + X2 + xf. 

We can now separate the displacement 6 = 6^ + 6j_ as follows. 

5|| =601 + 63 J3, 6± = 6iJi + 62 J2- (B.5) 

Now we are ready to study the first order term in the expansion, again 
using the representation of F{x) in terms of the exponential function. 

= Jdp f{p) p e^" ds e-^P" 6 e^^". (B.6) 

Since 6\\ commutes with x, the first part of this is simply F'{x) 6\\. For the 
part with 6j_ we use the formulas (lB.2p . where 6 is replaced by —spr. The 
integrals over s are trivial and we merely write sin{pr) and cos{pr) in terms 
of e^*^*" to get our final result. 

F{x + 6) = F{x) + F'{x)6ii + {F{x + ir) - F{x - ir)} — 6± + 

2ir 

{F{x + ir)+F{x-ir)-2F{x)} ^[J3,6^]+ 0{6^). (B.7) 
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It should be noted that the (5-related factors in Eq. (1B.7P can be written 
in the following way: 



[J3,Si] = ^[x,6], (B.8) 
S± = -\ [x,[x,6]], (B.9) 



6-6^. (B.IO) 



This means that we do not have to carry out the "rotation" that gave us Eq. 
(IB.4P explicitely; the talk about choosing a local coordinate system is merely 
rhetorical. 

I expect that this method can be extended to other Lie algebras, with the 
quantity S± subdivided into distinct portions according to the roots of the 
particular algebra. The system of Eqs. (IB. 91) . (]B.10|) would be adapted to 



make those separations, using the known values of the roots; and those root 
values would also appear in the final generalization of Eq. (]B.7p . 



Extension of this method to general matrix variables, over the complex 
numbers, is given in reference [S]. 

Appendix C - Fueter's Differential Equation 

The literature on Fueter's analysis of quaternionic functions points to 
a third order differential equation as his key result, extending the familiar 
Cauchy-Riemann (first order) equation for functions of a complex variable. 

Here I wish to present a rather simple derivation of that result, starting 
with the exponential function of our quaternionic variable. 

X = Xq + ixi + jx2 + kxs = Xo + u ■ X, u = {i, j, k), = x • x, (C.2) 

eP" = eP"«(cospr + u-x^^). (C.3) 

r 

Here is Fueter's first order differential operator. 

d . d . d , d d „ ,^ ^, 

OXq OXi 0X2 0X3 OXq 
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and we calculate its action on the exponential function and find the result 

□ = -2 eP^'^ (C.5) 
r 

This is not zero (as in the complex case, Eq. llC.ip ) but it is rather simple 
(and real). Moreover, we recognize this function as a solution of the four- 
dimensional Laplace equation. 



A4 = nn* = ^ + A3. (C.6) 

So here are two forms of Fueter's third order differential equation. 

A4 □ e^^ = 0, A4 □ = 0, (C.7) 

where the second result comes from expanding the first result in a power 
series in the parameter p. Thus any superposition (with real coefficients) 
of the powers or the exponential will satisfy this condition; and this is the 
basis for what they define as holomorphic functions of a quaternionic variable. 
They also exclude functions with terms such with arbitrary quaternion 

constants a - just as we have done in the present paper. 

From those differential equations (1C.7I) (Cauchy-Riemann-Fueter), some 
integral theorems follow. In the case of complex variables (IC.ll) . we get 
the familiar result that any integral around a closed path in the complex 
plane will be zero (with suitable analytic and single-valued behavior of the 
function f(z)). In the quaternionic case, the relevant integral is over a closed 
3-dimensional surface in the 4-dimensional space of the x^. 

Here is a surprise! Look at the result of our differential operator T) acting 
on the exponential function e^^, Eq. 04.51) . The coefficient of 5^ is the same 
function that we see on the right hand side of Eq. fIC.Sp . So, here is a new 
identification: 

dF{x) 1 



dx^ 2 
Appendix D - Caculating X> 



□F(x). (C.8) 



Here is a derivation of the last item in Eq. (18. 5p . which is made easy if we 
take a geometric perspective as we write the coordinate in four dimensions 
as X = Xq + rux- Start by writing 

V u^ = a5\\+ (3 5^, (D.l) 
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where a and f3 are to be determined. First, consider a displacement that has 
only S\\: this should not change at all, since such a displacement can only 
change xq and r. Thus, we see that a = 0. 

Next, consider a displacement that has only 6±: this should not change 
xq or r. So we write, 

V Ux = T>^ Ux = T>± {x — x^j/r = - V± x = - 6±. (D.2) 

Following the result for F^"^^ in Section 8, one may wonder whether the 
entire Taylor series might be written as, 

OO 1 

F(x + 5) = 5] - V'' F{x) = F{x). (D.3) 

fc=0 ^' 

This may be readily verified for the functions F{x) = x", starting with the 
equation, 

xe-^ = x + S. (D.4) 
For the exponential function, define 

Q{p) = e^^; (D.5) 

then calculate dQ/dp and use Eq. flD.4l) . 
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